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Abstract 



We analyze the LEP constraints from searches for charginos x an d neutralinos Xi> taking 
into account radiative corrections to the relations between their masses and the underlying 
Higgs-mixing and gaugino-mass parameters /i, and the trilinear mass parameter A t . 
Whilst radiative corrections do not alter the excluded domain in m x ± as a function of m x ± — 
m x , its mapping into the /i,mi/ 2 plane is altered. We update our previous lower limits on 
the mass of gaugino dark matter and on tan/3, the ratio of Higgs vacuum expectation values, 
in the light of the latest LEP data and these radiative corrections. We also discuss the 
viability of Higgsino dark matter, incorporating co-annihilation effects into the calculation 
of the Higgsino relic abundance. We find that Higgsino dark matter is viable for only a very 
limited range of /i and m^, which will be explored completely by upcoming LEP runs. 
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1 Introduction 



Considerable experimental effort is currently being devoted at LEP and elsewhere to the 
search for the charginos x an d neutralinos Xi j i = 1, .-,4 of the minimal supersymmetric 
extension of the Standard Model (MSSM) [Q, H, d, @|, and to understanding the constraints 
these and other particle searches impose on the parameter space of the MSSM || ||. One 
of the generic possibilities for MSSM phenomenology is that supersymmetry breaking 
is communicated to the observable sparticles via gravitational interactions, in which case 
the lightest supersymmetric particle is presumably the lightest neutralino Xi (hereafter x), 
which would be stable if R parity is conserved. The lightest neutralino is not directly visible 
to accelerator experiments in this type of /^-conserving supergravity scenario, and would be 
a good candidate @, || for the Cold Dark Matter advocated by many astrophysicists and 
theorists of cosmological structure formation [10], since it might well have 0.1 < Q x h 2 < 0.3, 
where Q x is the relic density of the lightest neutralino x i n units of the critical density, and 
h is the present Hubble expansion rate in units of 100 km/s/Mpc. Within this framework, 
which we adopt in this paper, detailed modelling of the relations between different searches 
is often necessary in order to close certain loopholes in the MSSM parameter space, in 
particular if one wishes to establish robust lower limits on the mass of the lightest neutralino 
X [0, [§• Such lower limits on m x are of particular relevance to the ongoing direct and 



indirect non-accelerator searches for cosmological relic neutralinos |L0 



The strength of the lower limit obtained depends || on various theoretical assumptions, 
such as the degree of universality among input supersymmetry-breaking parameters, and 
whether one requires that the relic density of neutralinos lie within the range favoured by 
astrophysical and cosmological considerations (K|. The successive upgrades of the LEP 2 
centre-of-mass energy |l|, [| |3|] have enabled this lower limit to be strengthened progressively, 
and it currently stands at m x > 40 GeV if the strongest versions of these theoretical and 
astrophysical assumptions are imposed M. 



Although radiative corrections to MSSM Higgs scalar masses |T|] have been included in 

the analysis of LEP data ||], and play an important role, radiative corrections to the masses 

of the charginos and neutralinos have not so far been included. The one-loop corrections 

to the relations between the physical masses and underlying MSSM parameters such as 

the gaugino masses M a : a = 1,2,3 and the Higgs- mixing parameter /i are by now well 

known [1^, [TJ], [TJ|] and non-negligible. [] In particular, they are known to modify considerably 

1 Radiative corrections to the x + X~ production cross section are also available H|, but appear less crucial 
at present, as we discuss below. 
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the tree-level difference between the lightest chargino and neutralino masses, m x ± — m x , 
which is one of the most important parameters controlling the efficiencies of experimental 
searches, especially in the "Higgsino" region of small \i and relatively large M a fT7| , [Hi . More 
generally, the radiative corrections to the chargino and neutralino masses cause shifts in the 
inferred exclusion domains in the /i, M 2 plane that are not negligible compared to those 
deduced from the improvements in the chargino mass limits obtained after successive LEP 2 
energy upgrades. These changes are also comparable to the differences between mass limits 
obtained under different theoretical assumptions. Therefore, radiative corrections to the 
chargino and neutralino mass relations could in principle play a role in the delicate interplay 
of the direct experimental constraints used to constrain m x [11, ||. Clearly the degree of 
effort put into the experimental searches |l], [|, |3| and their phenomenological interpretation 
merits the inclusion of these radiative corrections in the analysis of the LEP 2 data. 

We find that there are significant shifts in the allowed regions of the fi, Mi plane. The 
effects on m x ± are particularly important when it is analyzed in terms of M2, but not when 
it is analyzed as a function of AM = m x ± — m x . We use these radiative corrections in an 
update of our previous analysis of LEP constraints in the MSSM parameter space, and on m x 
and tan/3 in particular, assuming universal input soft supersymmetry-breaking masses and 
a plausible cosmo logical relic density. This update also includes the most recent preliminary 
limits on chargino and neutralino production from LEP running at 183 GeV |J1 [|, |J, 
as well as new limits on the lightest neutral Higgs mass rrih ]TP|, |2D|, [|]. Relaxing the 
universal scalar mass assumption, we focus on the portion of the MSSM parameter space 
for which the lightest neutralino is a Higgsino, where the limits are particularly sensitive to 



the radiative corrections [[17], [T]|. In this context, we evaluate the possibility of Higgsino 
dark matter, incorporating both co-annihilations and radiative corrections in conjunction 
with the bounds imposed by LEP. We find only a very limited range of parameters in which 
Higgsino dark matter is viable, which will be explored completely by forthcoming LEP runs. 



2 Implementation of Radiative Corrections in Chargino 
and Neutralino Codes for LEP 

As already mentioned, we work in the context of the MSSM in a supergravity framework 
where supersymmetry breaking is communicated to observable sparticles by gravity, so that 
the lightest neutralino \ is the lightest supersymmetric particle, and we also assume that 
R parity is conserved, so that x is stable. We further assume that the supersymmetry- 
breaking gaugino mass parameters are universal at the input supergravity scale: M a = mu2- 
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Similarly, we assume that the soft supersymmetry-breaking contributions to the squark and 
slepton masses rriq, m~ t are universal and equal to mo at the unification scale, and the trilinear 
parameters A are also taken to be universal. The MSSM can then be parameterized by 
m i/2, m o, A, (jl, the pseudoscalar Higgs mass m^, and tan (3. Here we project the experimental 
and cosmological bounds onto the fi, M 2 plane for fixed m , A, tua, and tan/?. 

We note that only the trilinear coupling A t , associated with the top-quark Yukawa cou- 
pling, is relevant for the phenomenology discussed here |, and that the possible range of val- 
ues of this parameter is constrained because of an infrared quasi-fixed point: A t ~ 2mi/2 ~ 
2rrig/3 [^l||, so that this parameter is inessential. This quasi-fixed point arises because, if A t 
is much larger than my 2, the leading-order running of A t is given by dA t /dt ~ 6XfA t /87r 2 [0, 
so that (for constant X t ) A t decreases as a power law with scale, with an exponent that is 
roughly 1/12 for X t ~ 1. Thus, if the input A t ^> my 2 , A t is reduced by an order of mag- 
nitude in its evolution from Mx to Mz- There are additional negative terms in the beta 
function for A t which are proportional to the gaugino masses, so that my 2 sets the scale 
for A t (Mz). Consequently, A t (M z ) 3> my 2 requires extremely large A t (M x )- However, 
this tends to drive the stop soft mass-squared parameters negative at Mz, so large values of 
At(Mz) are not allowed, at least for the values of tan/3 we consider. Therefore, we do not 
explore a range in values for A t , but concentrate on the fixed-point value. 

The input supersymmetry-breaking parameters are evolved down to the electroweak scale 
using the two-loop renormalization-group equations for the gaugino masses, so that the one- 
loop relation M a = (a a / olgut) x "^1/2 is violated at the 10% level. We also evolve the gauge 
and Yukawa couplings at the two-loop level, but other parameters are evolved with one-loop 
accuracy. At the tree level, we parameterize chargino and neutralino eigenstates in terms 
of M2, \i and tan/3. We take a top mass of 171 GeV: the only significant sensitivity of our 
results to m t comes in the Higgs mass constraints, which we discuss in Section 5. 

This theoretical framework has been adopted in many of the searches for supersymmetric 
particles at LEP, in particular for charginos and neutralinos. Until now, all the experimental 
analyses have used tree-level formulae for the masses and cross sections, in terms of the 
MSSM parameters introduced above. The one-loop calculations of the physical masses of 



the charginos and neutralinos |13|, [14], [TJ], and more recently the x + X production cross 



section [16], that have become available have been calculated in the DR prescription, which 



is the supersymmetric analogue of the conventional MS prescription. The one-loop mass 
formulae have now been incorporated in one of the codes available for LEP experimental 
analyses, and we discuss here a first exploration of their implications for sparticle mass limits. 



2 We do not consider very large values of tan/3 3> 10. 
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The radiative corrections to the cross section are less important for this purpose: they are 
generally < 15 % or so ||16|| , whereas the experimental upper limits on = a(e + e~ — > 
X + X~) are generally far below the theoretical predictions, except very close to the production 
threshold: Ecu ~ 2m x ± . Since varies very rapidly with m x ± in this region, the change 
the theoretical prediction for <r + _ has little effect on the chargino mass limit obtained. 

3 Exploration of Consequences of the Radiative Cor- 
rections 

Fig. |I] displays the effects of these one-loop radiative corrections in the fi, M2 plane. The thick 
lines are one-loop-corrected contours corresponding to fixed values of chargino and neutralino 
masses, and the thin lines are tree- level contours. The continuous lines are for m x ± =91 GeV, 
the dashed lines for m x ± = 86 GeV, corresponding respectively to the kinematic limits on 
m x ± from e + e~ — > x + X~ at the just completed 183 GeV run of LEP and its previous run 
around 172 GeV. The dash-dotted lines show the contour of m x = 40 GeV, corresponding 
to our previous lower limit on m x when several theoretical and cosmological constraints are 
applied 0. 

In general, each of the four neutralinos Xi is characterized by the four components 
a i;A;7i; and <5i of its eigenvector in the W , B, Hi, H 2 basis. In this notation, the gaugino 
fraction of a particular neutralino Xi is \J a i 2 + Pi 2 , and its Higgsino fraction is y 7i 2 + S 2 . In 
a later section of this paper, we discuss whether LEP limits still allow the lightest neutralino 
to be predominantly a Higgsino. With a view to this subsequent discussion, we recall here 
that the Z° coupling to a pair of neutralinos Xi and Xj is proportional to the Higgsino com- 
ponents of the neutralinos, namely 7i7j — 5i5j. Thus it is particularly easy for LEP to probe 
regions of the /1, M 2 plane where the lightest neutralino is mainly a Higgsino, by searching 
for production in association with a heavier neutralino state which is also mainly a Higgsino. 
Accordingly, we also plot in the relevant parts of Fig. [I] dotted lines representing the exper- 
imental limit of 0.18 pb on associated neutralino production e + e~ — > XiXj '■ i + j > 2. In 
the region of interest for this analysis, this is dominated by e + e~ — > XX2,3,4 and corresponds 
essentially to m x + m x ^ = 182 GeV, where x'h ls the lightest mainly-Higgsino state among 
the X2,3,4- 

The one-loop radiative corrections in Fig. Q are calculated with the following represen- 
tative values of the input MSSM parameters: tan/5 = 2, mo = 200 GeV, rriA = 1 TeV and 
At = 2m\/2- We see significant shifts between the pairs of corresponding thick and thin lines, 
making manifest the significance of the one-loop radiative corrections. Varying uiq between 
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Figure 1: The effects of one-loop radiative corrections in the //, M2 plane, calculated for 
tan(3 = 2, m = 200 GeV, uia = 1 TeV and A t = 2m\/2- The thick lines are one-loop- 
corrected contours corresponding to fixed values of chargino and neutralino masses, and the 
thin lines are tree-level contours. The continuous lines are for m x ± = 91 GeV, the dashed 
lines for m x ± = 86 GeV, the dash-dotted lines for the bound on associated neutralino produc- 
tion, and the dotted lines for m x = 40 GeV. 
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100 and 1000 GeV does not have a noticeable effect on the lines plotted, and they do not 
depend strongly on the precise value chosen for w,a, though they would be significantly dif- 
ferent for values of A t much larger than the fixed-point value that we favour. Interestingly, 
the differences between these pairs of lines are comparable to the splitting between the lines 
for m x ± = 86 and 91 GeV, indicating that the radiative corrections have phenomenologi- 
cal implications for the inferred constraints on the input MSSM parameters //, M 2 that are 
comparable to the upgrade in the LEP energy from 172 to 183 GeV. 

The LEP exclusions for charginos do not always extend to the kinematic limits. In 
particular, the lower limit on m x ± is sensitive to tuq: it affects cr + _ through the v exchange 
diagram, and when mo is such that m x ± — m„ < 3 GeV, chargino decays to a sneutrino 
and a low-energy lepton reduce the experimental efficiency to essentially zero. Moreover, 
as seen in Fig. |||, there is also a loss of sensitivity for fixed m at large M 2 , because of 
the reduction in detection efficiency that occurs when the difference in mass between the 
produced particle and its supersymmetric decay product, e.g., AM = m x ± — m x , is small. 
This implies that at large M 2 , beyond the top of Fig. [I], the experimental exclusion in the 



/i, M2 plane would contract toward ~ 45 GeV. It has been observed previously [1/, 18 



that AM may be altered significantly by the one-loop radiative corrections, altering the value 
of M 2 where this contraction would occur. For tan f3 not too large, the leading contribution 
to the change in AM, 5(AM), comes from heavy quark-squark loops, and for split stop 
masses it is approximately given by |T7], |TH[ 

3 „ max(m?, m~ ) 
S(AM) = —\l mt sin 29 t log X 2 . (!) 

0^7T TTLr 

ti 

where > m^ 2 . Note that this vanishes in the limit that there is no stop mixing (sin 26 t — > 
0), and that the sign of <5(AM) depends on the sign of 6 t . Both of these depend on the off- 
diagonal element of the stop mass-squared matrix, and hence on A t . Hence, the quasi-fixed 
point for A t is an important constraint on radiative corrections to AM. 

Incorporating these radiative corrections fully in the analysis of experimental data would 
require the use in experimental simulations of the full radiative-correction code P5|| . In the 
absence so far of such an analysis, we have used this code in conjunction with an ad hoc 
parameterization of the efficiency of the ALEPH chargino search as a function of AM to 
explore the possible significance of this effect. We see in Fig. §a that the exclusion contour in 
the M 2 ,m x ± plane is indeed affected significantly. However, as seen in Fig. ^p, this effect of 
the radiative corrections vanishes in the AM, m x ± plane, as expected from the small change 
in the chargino composition. This suggests that a useful rule of thumb may be to expect that 
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Figure 2: (a) The experimental limit on m x ± as function of M 2 and as a function of AM = 
m x ± — m x , for fixed itlq = 200 GeV and tan/3 = 2. The value of fi is determined by the 
combination m x ± and M 2 . The drop at large M 2 is due to the loss in experimental efficiency 
as x decays into x involve softer leptons. The dotted line comes from a tree-level analysis, 
and the solid line is obtained using an ad hoc parameterization of the experimental efficiency, 
in conjunction with the radiatively- corrected mass formulae. Note that the lines separate 
significantly when M 2 is large, (b) The experimental limit on m x ± is now plotted as a 
function of AM = m x ± — m x . Note that the tree-level and radiatively- corrected curves are 
almost coincident. 



the exclusion domain in this particular plane may be insensitive to the radiative corrections 



also in a complete experimental analysis using [23]. 



4 Updating Lower Limits on m x and tan/3 

We now review the lower limit on m x derived previously || using data from the LEP runs 
at energies up to 172 GeV, discussing two issues: the significance of the radiative corrections 
and the impact of the more recent data at 183 GeV. The purely experimental constraint from 
LEP 172 coming from chargino production and associated neutralino production yielded the 
limit m x > 14 GeV jl], with the minimum occurring around tan (3 = 3. Adding the theoretical 
assumption that all squark masses are universal at the unification scale, unsuccessful Higgs 
searches indirectly strengthened the bound on m x significantly in the region around tan (3 = 
2 |J. Furthermore, when scalar mass universality was extended to the Higgs sector and 
the values of fi and itla were fixed by requiring the consistency of electroweak symmetry 



breaking p2[ ) the bound on m x was improved for all tan/3. For tan/5 > 4, we found 
m x > 30 GeV, and the bound improved to m x > 40 GeV for tan (3 < 3 . 

We have also considered imposing the cosmological constraint that the neutralino relic 
density does not over-close the Universe. Specifically, the combination Q x h 2 is constrained 
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to be less than 0.3 by the requirement that the age of the Universe, to, be greater than about 
12 Gyr. These quantities are related through 



where the scaled Hubble parameter is h — ifo/(100km Mpc -1 s _1 ), and large relic densities 
may be incompatible with constraints on t - When the cosmo logical constraint was combined 
with the experimental constraints from Higgs searches and the theoretical assumption of 
universal Higgs masses, we found that values of tan (3 < 1.7 for \i < and < 1.4 for fi > 
were excluded. This is because the tree-level Higgs mass is small at low tan j3, so radiative 
corrections to must be enhanced by taking a large stop mass. On the other hand, large 
sfermion masses lead to large relic densities for a gaugino-like neutralino, as expected for 
universal scalar masses. Thus Q x h 2 < 0.3 is not compatible with a heavy enough Higgs boson. 
Finally, when we imposed the astrophysical desire that the cosmological relic neutralinos play 
a significant role in structure formation, which requires that Q x h 2 > 0.1, we found that the 
lower limit on the neutralino mass was strengthened at tan (5 > 3 to m x > 50 GeV. 

In the interests of simplicity, we re-discuss here the lower bound only in the case in 
which all theoretical and cosmological assumptions are made, in particular that the input 
scalar masses are universal and that the relic cosmological density lies in the range 0.1 < 
fl x h 2 < 0.3. The resulting bound on the LSP mass in this case was previously found to be 
m x > 40 GeV, with the minimum being achieved in a narrow range of 2.6 < tan (3 < 3.0, as 
seen in Fig. 1 of 0. In this case, the lower bound on m x was linked directly to the lower 
limit m x ± > 86 GeV from the LEP 172 GeV run and relied on the Higgs mass constraint. 
Cosmology did not play a direct role for this range of tan/3, though these constraints did 
play an important indirect role in closing off loopholes for special values of other MSSM 
parameters such as \x and mo where the limit on m x ± and hence m x could have been smaller, 
and provided stronger lower limits on m x for both lower and higher values of tan/3. 

We turn first to the impact on this analysis of the radiative corrections discussed in the 

previous section. We display in Fig. ^] the radiatively-corrected values of m x for m x ± = 

85 GeV, tan/3 = 2 and the representative value \i = —140 GeV, corresponding to m x = 

40 GeV at the tree level. The ranges plotted are generated by varying A t and mo over the 

ranges indicated in the two panels. We see that their effect is in general to reduce m x by ~ 1 

to 2 GeV. If we fix A t to the fixed-point value discussed earlier, we obtain the solid line 

shown in the second panel, and if we specify the value mo ~ 80 GeV that corresponded to 

3 We have verified that this shift is numerically insensitive to /i over a range larger than that of relevance 
to this analysis. 




(2) 
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the minimum value of m x in our previous analysis, we find a reduction in m x by ~ 1 GeV. 
We have investigated whether the roles of the LEP 172 GeV Higgs limits and of cosmology 
are altered significantly when radiative corrections are included, and found that this is not 
the case. We conclude that the sensitivity of the previous LEP 172 GeV lower limit on m x 
was around the 1 GeV level. 



M v = 40 GeV at tree level 




A 



1000 



m 



o 



Figure 3: Plots of the radiative correction to the lower limit on the neutralino mass, in the 
case that vn x ± = 85 GeV, tan(3 = 2 and fi = —140 GeV. The plots are obtained by varying 
-500 GeV < A t < 500 GeV and 50 GeV < m < 1000 GeV. The solid line in the second 
panel is for the fixed-point value A t = 2m 1 / 2 . We see that the limit on m x is reduced by 
about 1 GeV if we assume this value of A t 
tree-level lower limit on m. 



and take m = 80 GeV, corresponding to the 



We now turn to the implications of the LEP 183 GeV data set. Again, for simplicity we 
restrict our attention to the case of universal soft supersymmetry-breaking scalar masses, 
including Higgs masses, implement the cosmological density constraint and focus on the 
region tan/5 ~ 2 to 3, where we previously found the lower bound on m x under these as- 
sumptions. The LEP lower limit on the chargino mass is now m x ± > 91 GeV, corresponding 
to a strengthening of the lower limit on m x by ~ 3 GeV. Including the shift downwards 
by ~ 1 GeV due to the radiative corrections, this direct LEP lower limit corresponds to 
m x > 42 GeV. We have verified that the LEP Higgs searches and the cosmological density 
constraint continue to exclude loopholes in this limit. 

The best limit on the Higgs boson of the Standard Model is 88 GeV, which also applies 
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to rrih for tan/3 < 2 in the MSSM [[19], |20], [| [|]. The formulae we use to calculate rrih are 



believed to have an uncertainty of about 2 GeV |12|], so we allow MSSM parameter sets for 
which they yield rrih > 86 GeV at low tan/?. The LEP lower limit falls to about 78 GeV for 
tan/3 > 3, and we allow MSSM parameter sets for which m h > 76 GeV for such larger values 
of tan/?. These Higgs limits strengthen the lower limit on m x for tan/3 < 2.6 for ii < 0. 
Moreover, the strengthening of the LEP lower limit on m x ± , combined with cosmology, and 
including the two-loop running of the gaugino masses as well as the one-loop corrections to 
the chargino and neutralino masses, imply that m x > 50 GeV for tan/3 > 2.8. Only in the 
region 2.6 < tan/? < 2.8 for /i < is there a narrow cleavage where m x > 42 GeV is still 
possible. There is no corresponding cleavage for fi > 0, where the previous lower bound 
m x > 50 GeV for tan /3 > 2 remains essentially unchanged. 

As we discussed above, for sufficiently low tan/3, one cannot satisfy simultaneously the 
Higgs mass constraint and the requirement that to > 12 Gyr. With the strengthening of the 
lower limit on rrih from LEP 183, we now find that tan/3 > 2.0 for fi < and tan/? > 1.65 for 
fi > 0. Because the dependence of the Higgs mass on the stop masses is only logarithmic, 
these bounds on tan/3 are quite insensitive to variations in A t or to- These strengthened 
results make it very difficult to reconcile our assumptions of universal scalar masses and an 
interesting cosmological relic density with the low tan /3 infra-red fixed-point solution of the 
renormalizat ion-group equations favoured by some model builders |TI . 



5 Implications for Higgsino Dark Matter 

Although constrained supersymmetric models which impose universality on all the soft scalar 
masses, including the Higgs bosons, predict as a general feature that the lightest neutralino 



is a gaugino p4j as discussed above, it is well known that, if the universality assumption 
is relaxed, the lightest supersymmetric particle can be a Higgsino over roughly half the 
parameter space as viewed on the M 2 , [i plane ||, namely when < M 2 , as seen in Fig. f|, 
which displays contours of higgsino and gaugino purities, as well as radiatively-corrected 
m x contours. In this section we study whether it is still possible, in the light of radiative 
corrections and the latest LEP data, for the lightest neutralino to be a Higgsino, even if one 
relaxes the scalar-mass-universality assumption. 

For this purpose, our working definition of a Higgsino is that p 2 = 7, 2 + 5 { 2 > 0.81. In 
much of the Higgsino parameter space, it is possible to be even more definite. For small to 
moderate M 2 , the lightest supersymmetric particle is the Higgsino combination defined by 
S° = iJ°cos/? + sin/?, i.e., 7 = cos/? and 5 = sin/?, with m§ /i sin 2/3 JJ. Contours 
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Figure 4: Contours of neutralino purity: 99%, 97% and (for jj, < 0) 75%, and chargino 
and neutralino masses (solid lines). The long-dashed lines are contours of high bino purity, 
the dotted lines are contours of high photino purity, the dashed lines are contours of high 
H12 Higgsino purity, and the dash- dotted lines are contours of high So Higgsino purity. The 
radiative corrections to m x are calculated for tan/3 = 2 and itlq = 100 GeV. 



of 5*° purity are displayed in Fig. |] as dash-dotted lines. On the other hand, for large 
M 2 , the lightest neutralino is the state H\ 2 = -j%{H® ± -^2); i- e -> ^ = =^7 = il/v2 for 
sgn(/i) = ±1, with mi 



l h 12 



25 



Contours of H\ 2 purity are shown as dashed lines in 
Fig. |j. The cross-over point where the neutralino becomes more like if 12 than S° depends 
on jj, and tan /3, but is typically at M 2 of a few hundred GeV |25] , as seen in Fig. I The 
S° interacts much like a neutrino, i.e., it annihilates and scatters elastically via Z exchange, 
with a coupling reduced by cos 2 2/3. On the other hand, a "pure" H 12 state never annihilates 
or scatters through the Z, since the coupling 7 2 — 5 2 — > 0. Because of its similarity to a 
heavy neutrino, the S° may have an interesting relic density, with Q x h 2 determined mainly 
by its mass and tan (3. However, this same property makes it ideal for exploration by LEP. 
For example, for tan/3 = 2, a S° state with 99% purity exists only for \fi\ < 25 GeV, and 
from Figs. [I] and |] one can see that the region in the /x, M 2 plane for such a pure state is 
explored by LEP, as discussed below in more detail. In contrast, the H\ 2 has been very 
difficult to find or exclude. There are in principle two regions in the /i, M 2 plane where H\ 2 
can provide an interesting relic density. For large M 2 , they correspond to < My/ but with 
m x ± above the LEP limit, and > 1 TeV [p5|. The intermediate-mass Higgsino states are 



not of cosmo logical interest, because of their rapid annihilations to W and Z pairs. We have 
little to add concerning the very heavy H\ 2 states, but the lighter Higgsinos lie directly in 
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the region where the current LEP runs are eating away at the parameter plane, as we now 
discuss. 

In the region where the lightest neutralino is ifi 2 , this state is nearly degenerate with the 
associated chargino state, and in some cases with the second-lightest neutralino \2- There- 
fore, computation of the relic abundance must then take into account the co-annihilation 



with these degenerate partners |26], P7| . The relic abundance is determined in this case by 
a Boltzmann equation in which the annihilation cross-section is generalized to include co- 
annihilation terms: cx e ff = J2< J ij r i r j, where is the relative abundance of particle i, and we 
include all three nearly degenerate states. The factors are in fact exponentially sensitive 
to the mass differences between the states. These mass differences are in turn sensitive to 
the radiative corrections to the neutralino and chargino masses, with potentially significant 



implications for the relic density, as pointed out in [18 



As a first approximation, the relic densities calculated with and without radiative correc- 
tions are rather similar for similar values of AM. Deep in the Higgsino region at very large 
M 2 , the mass difference between the LSP and chargino is very small, and co-annihilations 
are very important in greatly reducing the relic density to uninteresting levels. In this limit, 
the small value of AM also makes this region inaccessible to LEP analysis. As one decreases 
M 2 , the mass difference increases and co-annihilation becomes less important, thus increas- 
ing the relic density and the Higgsino's cosmological significance, while at the same time 
making this region readily accessible to LEP particle searches. 

In Fig. |^ we display contours of constant chargino mass m x ± = 91 GeV, m x + rn x ' H = 
182 GeV, Higgs mass m^, and Higgsino purity, along with a contour of constant Q x h 2 = 0.1, 
for the choice tan/? = 2.0 which is of interest for the subsequent discussion. Since the 
universal scalar mass constraint must be relaxed for Higgsino dark matter, we have chosen 
rriA = rn Q = 1 TeV so as to maximize the Higgs mass and hence minimize the impact of this 
constraint. 

We see that the dashed lines in Fig. |5] representing the chargino and associated neutralino 
mass contours bound one away from small whilst the Higgs mass limit bounds one away 
from small M 2 . This is particularly restrictive at low tan/3, where the tree-level Higgs mass 
is small, and thus where radiative corrections to must be enhanced by taking large stop 
masses. The solid contour contains the region which leads to a significant neutralino relic 
density Q x h 2 > 0.1, and its limited range in M 2 is a result of co-annihilations. For larger 
values of M 2 , the neutralino is a purer Higgsino, and the masses of both the lightest chargino 
and the next-to-lightest neutralino approach the neutralino mass from above, enhancing the 
effect of co-annihilations that deplete the relic Higgsino abundance. For larger values of 
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Figure 5: Survey of experimental and cosmological constraints in the /i, M2 plane, focusing 
on Higgsino dark matter for tan(3 = 2 and (a) // < and (b) \x > 0. VFe pZoi t/je radiatively- 



corrected contours for m x ± =91 GeV, /or m x + 



182 GeV ; for selected values of nih 



and the Higgsino purity p, and for £l x h = 0.1. The shaded regions yield a Higgsino which 
satisfies the mass and relic density constraints described in the text. 



the relic density is suppressed by annihilations into W pairs. The hashed contours in Fig. |5|a 
represent Higgsino purities. Note the limited range of \i for which the mass and relic density 
constraints are satisfied. In Fig. [5|b we show the equivalent plot for \x > 0; in this case the 
Higgsino purities are lower, and the entire dark matter region falls below the purity cutoff. 

The combined effects of the above constraints, corresponding to the shaded regions of 
Fig. m are displayed for different values of tan/? in Fig. |6|. We have again taken mo = vtla = 
1 TeV to minimize the effects of the Higgs and chargino mass constraints. We find that there 
are no consistent Higgsino dark matter candidates for tan/? < 1.8 or > 2.5 for fi < 0, or for 
any value of tan/3 for /i > 0. The Higgs mass constraint cuts off the bottom of the allowed 
regions at low tan/5. When fi < it becomes a relevant constraint for tan/3 < 2.0 and is 
responsible for the complete disappearance of the allowed region when tan/3 < 1.8. Within 
the allowed regions displayed, the relic densities generally increase as \fi\ is increased, until 
the neutralino mass, whose minimum value here is ~ 71 GeV, becomes greater than mw, at 
which point the W + W~ annihilation channel opens, driving the relic Q x h 2 below 0.1. In any 
event, Q x h 2 is never greater than 0.12 anywhere in the allowed regions for fi < 0. 

If one considers neutralinos which have a lower Higgsino purity, i.e. are more mixed states, 
then additional allowed regions appear. For example, if the purity condition is relaxed to 
p 2 > 0.75, then for /i > there are regions of comparable area to those in Fig. |6| for which 



13 



300- 




ISO I 

-90 -80 -70 

Figure 6: The regions of the fi, Mi 'plane allowed by the constraints shown in the previous 
figure are shown for several different values of tanf3. There are no consistent choices of 
Higgsino parameters for tan (3 < 1.8 or > 2.5 for fi < 0, or for any value of tan j3 for fi > 0. 



Q x h 2 > 0.1, and these regions extend over a range in tan/? from 1.6 to about 25. Similarly, 
regions considerably smaller than those in Fig. ||] appear for fi < at tan /3 > 10 and persist 
up to tan/3 of about 45. At large tan/3, A t can vary somewhat from the quasi-fixed point, 
and as the Higgs mass is considerably larger than the current experimental bound, one has 
the freedom to play with mo to find regions of larger Q x h 2 . However, these regions remain 
very small, regardless of A t and rriQ. 

We have taken a large value of m , namely 1 TeV, for which the large stop masses produce 
large radiative corrections to m^. On the other hand, we have chosen a value of m t 2 
GeV below the current central value, and a larger top mass would increase the radiative 
corrections. However, this is important only for the low values of tan/3 discussed above, 
where the Higgs mass constraint is important. On the other hand, if tan/3 is too small, a 
larger top mass will lead to a Landau pole in the evolution of Aj. Due to this natural ceiling 
on the top mass as a function of tan /3, one cannot find any acceptable region for tan (3 < 1.7 
for fi < 0, even pushing the top Yukawa to its perturbative maximum. For /i > 0, we have 
already taken mt as large as possible, hence the allowed range of tan j3 cannot be extended 
to lower tan /3 in this way. 

Our allowed Higgsino regions shown in Figures || - ^| are somewhat smaller than those 



depicted in U18] , where the Higgsino region for negative fi extends up to M 2 ~ 400 GeV 
and down to ~ —60 GeV for Qh 2 > 0.1. An important difference in our analysis is the 
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Figure 7: Larger regions of the [i, M2 plane for tanf3 = 2 and (a) fi < 0, (b) fi > 0, with 
the radiatively- corrected contours for m x ± = 100 GeV and rrih = 100 GeV, three choices of 
the Higgsino purity, and the Q x h 2 = 0.1,0.3 contours for tuq = 100 GeV and (in FigQa) 
1000 GeV (dashed lines). We see that the allowed regions for gaugino dark matter are very 
much larger than those for Higgsino dark matter shown in Fig. [|. The small dark shaded 
region in Fig. ^b yields a stop with a negative mass squared. 



strengthening of the LEP bounds, not only on the chargino mass (to about 91 GeV), but on 
the associated production of Higgsinos (to m x + m x ' H > 182 GeV, as seen in Fig. ^). This 
alone cuts out about 2/3 of the Higgsino region of [ l8| . In addition, we have set the value 
of A t to its infrared quasi-fixed point: A t ~ 2mi/ 2 ~ 2nig/3 pi], for the reasons discussed 
above, whereas a larger value A t ~ 1 TeV was chosen in [[31|. Our choice of A t limits ab 
initio the size of the radiative corrections which could remove the degeneracy between the 
Higgsino and chargino states. Thus, in our analysis the co-annihilation of these particles 
remains important above M 2 > 250 GeV, in contrast to [pl| . 

The corrections to the chargino and neutralino masses are positive in the regime studied, 
and also increase with the sfermion masses, although the corrections are much smaller than 
for m/j. Thus the chargino and associated neutralino production bounds tend to be more 
restrictive for lower mo, unless one tunes so that it lies close to the chargino mass, in 
such a way that the production cross section is reduced by v exchange and the experimental 
efficiency for the chargino search is reduced. This might weaken the chargino constraint, but 
the associated neutralino constraint is actually strengthened by positive interference from 
v exchange. We also recall that the Q x h 2 contours themselves are quite insensitive to the 
choice of mo in the deep Higgsino region, so there is no great net change in the boundaries 
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of the allowed Higgsino regions in Figs. |5| and []. On the other hand, as we discuss below, in 
regions where the gaugino admixture is greater than in Figs. [5] and ^|, lowering mo decreases 
the relic density further as annihilations into fermions via sfermion exchange get enhanced. 
Lastly, as discussed above, the value of A t (mz) is quite insensitive to the choice of A t 

The areas A of the allowed Higgsino regions are never very large: calculating them 
using the logarithmic measure (d[i/ \i){dM2j M2) , we find that their maxima are reached for 
tan/3 ~ 2 for both /i < and /i > 0, with A max = 0.020, 0.006, respectively. Furthermore, the 
residual region for fi < could be explored by a search sensitive either to < 100 GeV or 
m x + m x / ff < 200 GeV, as anticipated if the LEP energy is eventually increased to 200 GeV, 
and the fi > region could be explored by a search for m x ± < 100 GeV 0. To put these 
numbers in perspective, one might wish to compare them with the logarithmic area of the 
full /i, M 2 plane for values of these parameters between 100 and 1000 GeV, namely A ~ 5. 
Alternatively, one might wish to compare with the logarithmic area of the parts of this plane 
where gaugino dark matter is still viable. This can be visualized in Fig. ^, where we display 
the regions where 0.1 < Q x h 2 < 0.3 for m = 100 GeV (solid contours) and m = 1000 GeV 
(dashed contours in Fig. ^a). Here we choose ttla = 3000 GeV to avoid the annihilation pole 
associated with the heavy Higgs, but the other parameter values are as in Fig. |5|. In the region 
with high gaugino purity and low Higgsino content, annihilation is predominantly through 
sfermion exchange into fermion pairs. For m = 1000 GeV, annihilation through this channel 
is effectively shut off, leading to the very small allowed region demarcated by the dashed 
contours in Fig. |S|a. For mo = 100 GeV, the cosmologically preferred region is significant. 
It is clear that the gaugino area is much larger than the remnant Higgsino area, indicating 
that this is a much more generic possibility. Indeed, the logarithmic area for gaugino dark 
matter is A > 2 for both [i positive and negative, some two orders of magnitude larger than 
the residual Higgsino regions. 

6 Conclusions 

The main purpose of this paper has been to include radiative corrections to chargino and 
neutralino masses for the first time in a phenomenological analysis of the LEP data. This 

4 We have checked what happens if we relax the gaugino mass unification constraint and take M\ = 
M2 |^| • We find that the allowed regions are somewhat smaller than in the unified case for /i < and 
disappear entirely for [i > 0. The latter results from a combination of larger neutralino masses - so that 
m x > rriw in a greater region, and closer ~ X mass degeneracy - leading to a higher co-annihilation rate. 

5 A fraction of the (very small) large tan (3 regions which open up for /1 < when the purity constraint is 
relaxed would not be probed by the LEP run at 200 GeV 
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analysis has been based on an implementation |23|] of the analytical results published pre- 
viously in one of the codes developed for LEP data analysis. This implementation makes 
it possible for radiative corrections to be taken into account in experimental analyses. In 
our first application of the radiatively-corrected code [23], we have found that radiative cor- 
rections generically reduce the upper limit on the lightest neutralino mass by about 1 GeV 
compared to what would be inferred from a tree-level analysis. 

We have also updated our previous analysis of constraints from LEP and cosmology 
to include preliminary results from the latest LEP run at 183 GeV. Assuming scalar-mass 
universality, the previous lower limit on m x is now increased to about 42 GeV, and the 
previous lower bounds on tan/3 are increased to about 2.0 for fi < and 1.65 for \x > 0, 
making life difficult for infra-red fixed-point models. 

We have also studied in detail the viability of Higgsino dark matter in the light of the 
latest LEP constraints, relaxing the scalar-mass-universality assumption. The radiative cor- 
rections play an essential role, as do the LEP constraints on charginos, associated neutralino 
production and Higgs bosons. We find that m x > 71 GeV in the remaining Higgsino region, 
that for jj, < there is only a very limited range of tan/3 over which Higgsino dark matter 
is viable at all, and that the area of the the /i, Mi plane within which it is viable is always 
very restricted. This is particularly evident when the Higgsino region is compared with the 
region of the fi, M 2 plane over which gaugino dark matter is viable. 

Finally, we review briefly the prospectives for further exploration of neutralino dark 
matter with higher-energy runs of LEP. We have already observed that prospective im- 
provements in the sensitivity to MSSM Higgs bosons, charginos and other supersymmetric 
particles should enable the complete investigation of neutralino dark matter up to ~ 50 GeV, 
if universal input scalar masses and an interesting cosmological relic density are assumed. 
Additionally, the analysis of this paper indicates that the future higher-energy LEP runs 
should also provide closure on the possibility of Higgsino dark matter, by telling us whether 
the lighter chargino or the lightest MSSM Higgs boson weighs less than 100 GeV. 
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